Cooling and squeezing via quadratic optomechanical coupling 
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We explore the physics of optomechanical systems in which an optical cavity mode is coupled parametrically 
to the square of the position of a mechanical oscillator. We derive an effective master equation describing two- 
phonon cooling of the mechanical oscillator. We show that for high temperatures and weak coupling, the steady- 
state phonon number distribution is non-thermal (Gaussian) and that even for strong cooling the mean phonon 
number remains finite. Moreover, we demonstrate how to achieve mechanical squeezing by driving the cavity 
with two beams. Finally, we calculate the optical output and squeezing spectra. Implications for optomechanics 
experiments with the membrane-in-the-middle geometry or ultracold atoms in optical resonators are discussed. 
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Introduction. In optomechanical systems optical and me- 
chanical degrees of freedom are coupled via radiation pres- 
sure, optical gradient, or photothermal forces. While work in 
this area was originally motivated by the goal of building sen- 
sitive detectors for gravitational waves the field has be- 
come an active area of research in its own right. Its main goal 
is to investigate quantum coherence in macroscopic solid-state 
devices both for quantum information purposes and gaining 
new insights into the quantum-to-classical transition (3j|4l]. 

In most optomechanical experiments an optical cavity 
mode is parametrically coupled to the position of a mechanical 
oscillator. Consequently, many properties of this setup have 
been discussed, including red-sideband laser cooling in the 
resolved-sideband limit ||5l0], normal-mode splitting 10,181], 
optical squeezing llHGH], backaction-evading measurements 
IUll - [l3ll . mechanical squeezing using either feedback IU2l[l4ll . 
squeezed light |[T5ll or modulation of input power llrjll . and en- 
tanglement between light and a mechanical oscillator [IT]. 

However, some optomechanical systems feature a quadratic 
optomechanical interaction, i.e. an optical cavity mode is cou- 
pled parametrically to the square of the position of a mechan- 
ical oscillator. One example is the membrane-in-the-middle 
geometry, in which the membrane is placed at a node or antin- 
ode of the cavity field II 1 8142011 . A second system which can 
realize quadratic optomechanical coupling is a cloud of ultra- 
cold atoms loaded into an optical cavity, where the cloud's 
center-of-mass coordinate serves as the mechanical degree of 
freedom I2H l22ll . To date, the theoretical literature has fo- 
cused on usin g the qu adratic coupling to detect phonon Fock 
states B18[|l9ll23l425ll . but otherwise the possible uses of this 
form of optomechanical coupling are largely unstudied. 

In this paper we explore three features of quadratic optome- 
chanical coupling: two-phonon cooling of the mechanical os- 
cillator, squeezing of the mechanical oscillator, and squeez- 
ing of the optical output field. Using Fermi's Golden rule we 
first write down an effective master equation for the mechan- 
ical oscillator. In the classical limit of large phonon number, 
two-phonon cooling processes change the steady-state num- 
ber distribution from exponential to Gaussian, a consequence 
of the nonlinear damping. In the quantum limit we find that 
ground-state cooling is not possible since two-phonon cool- 
ing processes preserve the phonon-number parity. We then 
demonstrate that the model maps onto a degenerate paramet- 



ric oscillator if the cavity is driven by two laser beams whose 
frequencies are detuned to either side of the cavity resonance 
by an amount equal to the mechanical frequency. This opens 
up the possibility of mechanical squeezing. Finally, we cal- 
culate the optical output spectrum and find that this system is 
capable of producing optical squeezing. 

Hamiltonian. We start from the Hamiltonian (with H = 1) 

H = (wr + gx 2 ) (a) a - {a) a)) + u> M &b + H y + H K (1) 

where ujr is the cavity resonance frequency, g the quadratic 
optomechanical coupling, and x — xzPF{b + b*) the posi- 
tion of the mechanical oscillator with zero-point fluctuations 
xzpf — (Smwjvf) -1 / 2 , frequency lum and mass m. a and b 
are annihilation operators obeying bosonic commutation re- 
lations, and H K describe the coupling to the mechan- 
ical and optical baths and the optical drive. We have sub- 
tracted the steady-state mean photon number (d^a) which 
renormalizes the frequency of the mechanical oscillator. The 
Hamiltonian (fl~|i is relevant to systems with membrane-in-the- 
middle geometry [18], to ultracold atoms in optical resonators 
lHHim, and double-microdisk whispering-gallery mode res- 
onators rt33tl when the first derivative of the cavity disper- 
sion relation w cav (a;) vanishes, i.e. uj' c&v (xq) = 0, so that 
g = w" av (a;o)/2 is the leading order of the optomechanical 
coupling. 

Expressing a = e~ lulLt (a + d) with the laser frequency wi, 
choosing a real, and neglecting d?d with respect to d(d' + d), 
we obtain the following quantum master equation 

Q=-i Hs,Q + nV[d]Q + 7(1 + n th )T>[b]g + 7 7i th 2?[6 t ]g 

(2) 

with the system Hamiltonian 

H s = -ASd + u M Pb + g(b + tf) 2 (S + d) (3) 

where A — lul — uir is the detuning, g = gax^pp the cou- 
pling, k and 7 the cavity and mechanical damping rates, and 
7i t h the thermal phonon number. We assume that the optical 
bath is at zero temperature. 2?[6]p = ogo^ — (o^og + gd^d)/2 
denotes the standard dissipator in Lindblad form. 

Current experiments with the membrane-in-the-middle ge- 
ometry work in the resolved-sideband limit at large thermal 
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FIG. 1: Steady-state mean phonon (wo) (left) and photon number 
(cftd) (right) as a function of the coupling strength g/K and thermal 
coupling j/k for thermal phonon number nth = 1, obtained from 
the numerical solution of the full quantum master equation l[2}. 

phonon number, i.e. k = 10 5 Hz, u>m — 10 6 Hz, 7 = 0.1Hz, 
and nth = 10 7 at T = 300K, while the optomechanical cou- 
pling is small compared to the cavity linewidth g/n = 10~ D 
IU8M20I1 . Experiments with ultracold atoms in optical res- 
onators have recently also realized quadratic optomechanical 
coupling, albeit with much larger coupling g/n s» 1 at small 
thermal phonon number n t h- However, they operate in the 
single-photon regime and outside the resolved-sideband limit, 
so that our results do not directly apply 112 ll 12211 . 

Numerical simulations of the full quantum master equa- 
tion. If the cavity is driven on the red two-phonon resonance, 
i.e. A = — 2u>m> m the good-cavity limit, i.e. k -C we 
expect two-phonon cooling processes to be important. Con- 
centrating on the resonant terms, i.e. Hs = giwwd + h.c), 
we solve the full quantum master equation (O numerically for 
small thermal phonon numbers n t h- As the coupling strength 
g increases, we find that (i) for jn t h <C k the mechanical os- 
cillator is cooled due to the coupling to the zero-temperature 
bath of the optical field l26ll . (ii) for 771th ^> K the optical 
field is heated by the coupling to the mechanical oscillator, 
and (iii) for 7n t h ~ k both effects are important and the den- 
sity matrix has non-zero off-diagonal elements so that e.g. the 
correlator | (d* bb) | is non-zero. These features are generic, but 
the details depend on the thermal phonon number n t h- As an 
example we plot in Fig.[T]the steady-state mean phonon num- 
ber (wb) and the mean number of photons due to the coupling 
to the membrane (d'd) as a function of the coupling strength 
g/n and thermal coupling j/k for n t h = 1. 

Effective master equation describing two-phonon cooling. 
For 7nth "C k and weak coupling, i.e. j < c, we can em- 
ploy a quantum noise approach M26I1 . i.e. we calculate two- 
phonon cooling and amplification rates using Fermi's Golden 
rule. Concentrating on the diagonal terms of the density ma- 
trix g nn — P n , we write down a set of rate equations 

P n = -7 (ntb(n + 1) + («th + 1)") Pn 
+ jn th nP n -i +7(n th + l)(n+ l)P„+i 
- (I>(n - 1) + r t (n + 2)(n + 1)) P n 
+ r 4 (n + 2)(n + l)P„ +2 + r t n(n - 1)P„_ 2 . (4) 

The terms in the first two lines are due the coupling of the me- 
chanical oscillator to its thermal bath with rate 7 and thermal 
phonon number n t h- The terms in the last two lines are due 



to two-phonon processes whose rates are governed by = 
g 2 x^ PF S nn (±2uj M ) where S nn (uj) = n\a\ 2 \x R (u)\ 2 is the 
photon number spectral density and = [k/2 — i(oj + 

A)] -1 the cavity response function. We note that nonlinear 
damping of an oscillator has been studied in Ref. j27il28ll . 

The infinite set of rate equations (01 can be replaced 
by a single differential equation for the generating function 
F(z, t) = X)^Lo Pn(t)z n . For vanishing two-phonon ampli- 
fication r<|- = the steady-state equation can be solved ex- 
actly in terms of the confluent hypergeometric function 12911 . 

We point out that it is an advantage of quadratic cooling that 
it enables cooling when the membrane or atoms are placed at 
a node of the cavity field. In this situation the system is most 
insensitive to absorption from the membrane or atoms as well 
as to bistability from radiation pressure. 

Two-phonon cooling in the classical limit. For large ther- 
mal phonon number nth 3> 1 and weak two-phonon cooling 
7n t h ^> T^ we can replace the quantum operators b and d 
in their Heisenberg equations of motion by complex ampli- 
tudes j3 = (b) and a = (d) and obtain two coupled clas- 
sical Langevin equations j3 — —j/3/2 — 2ig(3*a + £ and 
a = —na/2 — ig/3 2 , where the thermal noise is characterized 
by (£(*)£*(*')) = JTithS(t—tf). Adiabatically eliminating the 
optical field we obtain an equation of motion with a nonlinear 
damping term /3 = -7/3/2 - 4g 2 \/3\ 2 (3/ k + £. Solving the 
corresponding Fokker-Planck equation we obtain the phonon 
number distribution 
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The distribution changes from an exponential for 7/r^ 3> n t h 
to a Gaussian for 7/T1 <C n t h- In the latter limit the mean 
phonon number is given by (n) = y '"/n t hn / 'irg 2 '. The ex- 
pression in Eq. (0 agrees with the high-temperature limit of 
the exact solution to the rate equations II29I1 . We conclude that 
the change in the steady-state phonon number distribution is a 
purely classical effect due to nonlinear damping. 

The quadratic coupling in current membrane-in-the-middle 
experiments is very small. Nonetheless, it leads to sizable ef- 
fects if the thermal phonon number n%h is large. In Fig. [2] 
(left) we plot the steady-state phonon number distribution for 
a thermal phonon number n t h = 10 7 , corresponding to a tem- 
perature of T — 300K, both in the absence — and pres- 
ence r^/7 = 4 ■ 10~ 7 of two-phonon cooling. We find the 
Planck distribution with mean nth = 10 7 becomes a nearly- 
Gaussian distribution with mean nth = 2.8 • 10 6 . The qualita- 
tive change in the phonon number distribution for membrane- 
in-the-middle experiments is the main result of this paper. 

Two-phonon cooling in the quantum limit. Solving the rate 
equations © in the limit of strong optical damping and/or 
small thermal heating rate 771th. i.e. Tj, 3> 771th, we ob- 
tain the phonon number distribution Pq/P\ = 3 + 1/nth 
with all other P n = and the minimal mean phonon number 
n M = 1/(4+1 /n t h). The fact that strong cooling leaves both 
ground and first excited state occupied is a consequence of the 
fact that two-phonon cooling processes preserve the phonon- 
number parity. 
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FIG. 2: (Left) Classical limit: Semi-logarithmic plot of the steady- 
state phonon number distribution P„ for thermal phonon number 
Thh = 10 r . Planck distribution for Tj_ = (dashed) and nearly- 
Gaussian distribution for r^/7 = 4 ■ 10 -7 (solid). (Right) Quantum 
limit: Mean phonon number (6 6) as a function of coupling qjn 
obtained from the quantum master equation (blue dots) and the clas- 
sical rate equations without (red solid) and with strong-coupling cor- 
rection (green dashed). The black dashed line indicates the minimal 
phonon number n^f f° r n th = 5. 



In Fig. |2] (right) we plot the steady-state phonon number 
(wb) as a function of coupling g obtained from the full quan- 
tum master equation (O and the rate equations (|4j. The two 
results show excellent agreement at weak coupling g <C n. 
As 4c? 2 /k becomes comparable to k, the Fermi's Golden rule 
expression for two-phonon cooling breaks down and the pre- 
dictions of the full master equation (O and the weak-coupling 
rate equations will in general be different. In the limit 
where only the Fock states with phonon number n < 3 are 
important, we adiabatically eliminate the off-diagonal terms 
in the quantum master equation (0 and find that the strong- 
coupling two-phonon cooling rates are given by 



2 n(n — 1)k 



+ 4g 2 n(n — 1) 



(6) 



In the weak-coupling limit g <C n this expression simplifies 
to our previous result Ag 2 n(n — 1)/k. In the limit of strong 
coupling g ^> k, the two-phonon cooling rate n _ 2 remains 
finite and cannot exceed the cavity damping rate k. This leads 
to a minimum phonon number which is larger than the one 
predicted from the weak-coupling theory. In Fig. [2] (right) we 
plot the steady-state mean phonon number obtained from the 
modified rate equations and find excellent agreement with the 
exact solution to the full quantum master equation ©. 

We note that for quadratic coupling in the good-cavity limit, 
the phonon number distribution P n can be measured by mon- 
itoring the phase shift of the reflected light II la. 1 19112311 . 

Mechanical amplification and squeezing. Driving the cav- 
ity at both usr ± u>m with equal strength, the classical part of 
the cavity field oscillates in time a = Acosu>Mt. In the case 
of linear optomechanical coupling this enables a backaction- 
evading measurement of one quadrature of the mechanical os- 
cillator UH-Q"!], but does not itself produce squeezing apart 
from the one which is conditioned on the measurement out- 
come. When we instead consider the same drive applied to a 
system with quadratic coupling, moving to an interaction pic- 
ture with respect to Hq — LORC^a + LJM^b and keeping only 
non-rotating terms of the quadratic coupling, we obtain the 




FIG. 3: Optimal squeezing spectrum of the cavity output field 
SopJ^ixj) asafunction of detuning A. The parameters are ui m /n = 5, 
g/n — 0.5, j/k = 0.1 and n t h = 0. 



standard Hamiltonian of the degenerate parametric oscillator, 
i.e. H DP0 = § (P + (b f ) 2 ) with X = gxl PF A/2. Thus this 
setup could be used to amplify small mechanical signals l26ll . 

Solving the linear quantum Langevin equations we see that 
the steady-state fluctuations in the quadrature X — (be lUMt + 
H.c.)/ \pl of the mechanical oscillator are squeezed below the 
thermal level, (X 2 ) = "^^x/y > depending on the ratio x/l- 
At threshold \ = 7/2, i.e. before the parametric oscillator 
becomes unstable, the squeezing reaches the theoretical limit 
for an internal mode of 3dB 03OI1 . Although the parameters of 
current membrane-in-the-middle setups show small coupling 
g <§; k and large thermal phonon number nth 3> 1, we empha- 
size that the ratio xH which is of importance here can still be 
comparable to unity and lead to significant noise squashing. 

The mechanical squeezing can be detected by coupling the 
position of the mechanical oscillator parametrically to a sec- 
ond optical mode. This scenario has been studied in Ref. Hill . 

Output spectrum. Let us now return to the full model OJ. 
Up to second order in the coupling g, the cavity output spec- 
trum which is defined as S^M = / dt e tut (dt ut (t)d out (0)) 
with d out = d in + y^nd and (d in (t)d} n (t)) =5(t-t') is given 
by Sgf (w) = AKg 2 \ XR [-u]\ 2 S x * x z{u) where 
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7(nth + 1) 2 
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l n th 



2jn t h(n t h + 1) 



7 2 + (w + 2lu m ) 2 7 2 +w 2 



(7) 



We see the output spectrum Sj^fa) h as sidebands at to = 
±2uj m and w = as expected for quadratic coupling. 

Optical squeezing spectrum. One application of optome- 
chanical devices which has been widely advocated |0, [Toll 
is to use the nonlinear coupling between light and mirror to 
squeeze the incoming coherent light beam, i.e. reduce one of 
its quadratures below the shot-noise level at certain frequen- 
cies. The quantity which characterizes this noise reduction is 
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the optical squeezing spectrum Sg ut (uj) given by 13111 

/>oo 

S% ut (uj) = 1+1 dte lu)t {: X° ut (t) , X° ut (0) :) 

Xg(t),X (O) :> 



= 1 + K dt e lut T 



where X° ut = (dl ut e* +H.c.)/2, (A, B) = (AB)-(A)(B), 
the colons indicate normal ordering and T is the time-ordering 
operator ll3lll . The former expression is useful for calculations 
in input-output theory and the latter for master equation sim- 
ulations evoking the quantum regression theorem. 

In Fig. [3] we plot the optimal squeezing spectrum of the 
cavity output field S°^\(ui) = mine Sg ut (to) as a function 
of detuning A. In the good-cavity limit at zero temperature 
we find squeezing at oj = and u> = 2lom for A = and 
A = — 2lum, respectively. We choose parameters u m / 1 k = 5, 
g/n = 0.5, j/k = 0.1 and n t h = 0, possibly relevant to 
future experiments with ultracold atoms in optical resonators. 

To gain insight beyond numerics we obtain the squeezing 
spectrum perturbatively up to second order in the coupling g 

s$M = 1 + s^TH + s° d f(~u) 



Xr{u)Xr{~u)S X 2 X 2 (uj) 



1 



«th + - ) X*r{u)Xr(-u) (kxr(u) - 1) 



7 + i(w — 2um) 7 + i(w + 2wm) 



As in the case of linear optomechanical coupling UQli, we 
find two regimes of squeezing: for small detuning and small 
frequencies as well as for detuning and frequencies close to 
twice the mechanical frequency. For the same parameters and 
at zero temperature the amount of squeezing is comparable to 

(8) the linear optomechanical systems. However, thermal fluctua- 
tions are more destructive for squeezing effects with quadratic 
optomechanical coupling due to the quadratic scaling with the 
thermal phonon number nth and due to the fact that the spec- 
trum S x i x i (oj) has weight at zero frequency. 

Conclusions. We have explored the physics of nonlinear 
optomechanical systems where an optical cavity mode cou- 
ples quadratically rather than linearly to the position of a 
mechanical oscillator. For optomechanical experiments with 
membrane-in-the-middle geometry, we predict a qualitative 
change in the phonon number distribution and mechanical 
noise squashing. For future experiments with ultracold atoms 
in optical resonators, we found the quantum limit of two- 
phonon cooling as well as mechanical and optical squeezing. 
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